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The Dirac electrons of graphene, an intrinsic zero gap semiconductor, uniquely carry spin and 
pseudospin that give rise to many fascinating electronic and transport properties. While isolated 
zigzag graphene nanoribbons are antiferromagnetic semiconductors, we show by means of first- 
principles and tight-binding calculations that zigzag graphene nanoribbons supported on graphene 
are half-metallic as a result of spin- and pseudospin-symmetry breaking. In particular, half-metallic 
Dirac cones are formed at K (K’) near the Fermi level. The present results demonstrate that 
the unique combination of spin and pseudospin in zigzag graphene nanoribbons may be used to 
manipulate the electronic properties of graphene, and may have practical implications for potential 
graphene-based nanoelectronic applications. 

PACS numbers: 71.20.-b,73.20.-r,73.22.Pr,75.75.-c 


I. INTRODUCTION 

The discovery of graphene has inspired intense inter¬ 
est in graphene-based nanostructures and their deri va- 
tives such as van der Waals (vdW) heterostructures .1^^ 
One fascinating feature of graphene is the combina¬ 
tion of real spin and pseudospin which gives rise to 
rich physics in this unique two dimension system, e.g., 
unusual quantum Hall effect .1^ Breaking these symme¬ 
tries could lead to interesting physical phenomena, 
such as gap opening due to pseudospin symmetry break¬ 
ing in supported graphene,^ additional Hall plateau in¬ 
duced by SU(4) spin-pseudospin symmetry breaking due 
to electron-electron interaction!^ and valley polarization 
caused by inversion symmetry breaking.^ 

Graphene can be patterned into a variety of one di¬ 
mensional nanoribbons, of which electronic properties 
are strongly dependent on the ribbon width and edge 
patternsAmong them, zigzag graphene nanorib¬ 
bons (ZGNRs) have attracted particular attention ow¬ 
ing to their uniqne edge magnetism,which can be 
strongly suppressed by metal substrates.l^^M^ Recently, 
H-terminated ZGNRs on graphene substrate were ob¬ 
tained by cutting the top layer of graphene bilayer us¬ 
ing hydrogen etching method.l^ Because of the van der 
Waals interaction between the GNRs and the substrate, 
this type of system provides an ideal platform for study¬ 
ing electronic properties of GNRs. For instance, edge 
states were explicitly observed in STM/STS experiments 
and a critical width of about 3 nm was revealed for the 
onset of electron-electron correlation between the edges 
of anti-ferromagnetic (AFM) ZGNRs.^^ 

The edge states of ZGNRs also carry unique spin and 
pseudospin information: spins are ferromagnetically cou¬ 
pled along each edge but antiferromagnetically coupled 
between edges, i.e., along each edge the spin polarized 
state is spatially located on only one sublattice (pseu¬ 
dospin). Here we demonstrate, based on density func¬ 
tional theory (DFT) and tight binding (TB) calculations, 
that the unique combination of spin and pseudospin in 
ZGNRs can also be used to manipulate properties of 



FIG. 1. Top view of AB stacking hydrogen-terminated AFM- 
ZGNR/graphene. Green and brown balls represent carbon 
atoms in the nanoribbon and graphene, respectively. Red 
arrows denote spins on the edge atoms. The black box rep¬ 
resents the unit cell for calculations. A ZGNR with n zigzag 
chains is denoted by n-ZGNR. The nanoribbon here is 8- 
ZGNR. 


graphene. A half-metallic state is surprisingly obtained 
in AB (Bernal) stacking of ZGNR on graphene as a re¬ 
sult of spin- and pseudospin-symmetry breaking caused 
by the interaction between the GNR and graphene sub¬ 
strate. Remarkably, half-metallic Dirac cones are formed 
at K (K’) near the Fermi level. 


II. COMPUTATIONAL DETAILS 


To study the electronic and magnetic properties 
of AFM ZGNRs supported by graphene, we carried 
out spin-polarized first-principles electronic-structure 
calculat ions using the Vienna Ab initio Simulation 
Package! ^^ * ^^ ! The exchange-correlation functional was 
parametrized in terms of the local density approxima¬ 
tion according to Ceperley and Alder,!^ and pseudopo¬ 
tentials were const ructed by the projector augmented 
wave method.^^lllll xhe one-dimensional Brillouin zone 
was sampled by a Ix24x 1 Monkhorst-Pack mesh for the 
self-consistent calculations. 
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FIG. 2. DFT-derived electronic bands for AA stacking of 8- 
ZGNR/graphene. The Fermi level is represented by a dashed 
line. 


III. RESULTS 
A. DFT calculations 

Figure [l] schematically illustrates the structure of Id- 
terminated AFM-ZGNR/graphene with AB stacking. 
The ribbon is separated from its periodic images by ~12 
Avacuum regions. In both the DFT and the TB calcu¬ 
lations the spin up (down) density corresponds to the 
majority spin of the left (right) edge of the nanoribbon, 
which is over the top (hollow) sites of the graphene sub¬ 
strate. Our DFT calculations indicate that edge mag¬ 
netism of AFM-ZGNRs is preserved in the presence of 
graphene, with a magnetic moment ~0.13 fis per edge 
atom for both the freestanding and supported ZGNRs. 
Although the magnitudes for the two edge are slightly dif¬ 
ferent for the supported system, the difference is small; 
for example, for 32-ZGNR/graphene, the two edges have 
moments of 0.128 and 0.136 fiB, respectively. (The mag¬ 
netic moments on the edge atoms obtained in the present 
study are much smaller than the one (0.43 fis) reported 
by Ref. [TTl but consistent with Ref. [ini However, we 
note that larger moments (~0.4 /i^) can be obtained if 
the core correction to the density is neglected.) 

To see the effects of stacking on the electronic structure 
of ZGNRs on graphene, calculations were also performed 
for AA stacking 8-ZGNR/graphene. As shown in Fig. ^ 
AA stacking induces a small gap of about 0.03 eV at the 
K point. The two spins are degenerate since A A stacking 
maintains the inversion symmetry of the nanoribbon. 

Figures ir and d display the electronic bands for 8- 
and 32-ZGNR/graphene, respectively. Compared to AA 
stacking, a spin splitting is obtained for AB stacking 
ZGNR/graphene. The gap opening that occurs at the 
projection of the K point (27r/3a) of graphene is spin- 
dependent: it is much larger for spin-up - the major¬ 
ity spin of the left edge (on top sites) - than that for 
spin down, the majority spin of the right edge (in hol¬ 
low sites). In particular, for spin down the gap between 
the valence band and the conduction band is negligibly 


small at K. The trend in the gap opening is in fact associ¬ 
ated with the location of the edge atom (top sites versus 
hollow sites). The Fermi level crosses bands of only one 
spin state, making the whole system half-metallic as a re¬ 
sult of the nanoribbon-graphene interaction that induces 
small magnetization in the graphene and a net magnetic 
moment to the whole system. 

Figures [^, c show the electronic bands for 8- 
ZGNR/graphene weighted by the localization on the 
graphene substrate and the nanoribbon, respectively; 
Figs. n and f are the corresponding plots for 32- 
ZGNR/graphene. One can see that the bands for 8- 
ZGNR/graphene at K near the Fermi level are basi¬ 
cally graphene bands. For 32-ZGNR/graphene, both 
constituents have considerable contributions, as expected 
since the bands evolve toward the properties of graphene 
bilayer as the ribbon size increases. Band analysis fur¬ 
ther reveals that the four bands (including spin) at K are 
contributed by the hollow sites only, consistent with the 
limit of graphene bilayer. However, instead of the four¬ 
fold degeneracy at the K point in a graphene bilayer, for 
AFM-ZGNRs/graphene the four bands split as a result of 
spin and pseudospin symmetry breaking. From Figs. n 
and f one can see a shift in the edge states near the X 
point, that is, the bands for the left edge (top sites) are 
shifted upward relative to those for the right edge (hollow 
sites). The underlying physics is that the AB stacking of 
AFM-ZGNRs and graphene gives rise to different elec¬ 
trostatic potentials for different edges: Because the edge 
states are mainly localized to one sublattice, especially 
near the X point, they are either at the top or hollow 
sites. The AB stacking raises the electrostatic potential 
of the top sites more than the hollow sites, resulting in a 
relative shift in the bands for the two edges. This mech¬ 
anism is similar to the case of freestanding AFM ZGNRs 
under external electric field.lSl Dramatic changes can be 
seen as the k point approaches K, where the shift of the 
conduction bands of the two spins is opposite to the case 
near the X point. This behavior is related to the inter¬ 
action between states of AFM ZGNRs and graphene in¬ 
volving special spin- and pseudospin-symmetry breaking 
(see discussions below). 

To further explore the properties of the bands at 
K (K’), Fig. 1^ depicts the ^-projected bands around 
K for both 8-ZGNR/graphene and 32-ZGNR/graphene, 
obtained by projecting the wave functions of AFM- 
ZGNR/graphene in the r ectangular supercell onto the 
1x1 cell of graphen^sSEI] ^ pjg Half-metallic 

Dirac states about K exist for all the studied sys¬ 
tems. However, unlike for ideal graphene, the Dirac 
cone for ZGNRs/graphene is anisotropic, i.e., the lin¬ 
ear dispersions along F-K and its perpendicular direc¬ 
tion are different, leading to anisotropic Fermi velocities 
in ZGNRs/graphene due to the one-dimensional nature 
of the nanoribbon. Moreover, the linear dispersions for 
32-ZGNR/graphene experience a significant renormaliza¬ 
tion, as expected since the bands about K approach those 
of a graphene bilayer. (Galculations were also performed 
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FIG. 3. DFT bands for AB stacked AFM-ZGNR/graphene for (a) for 8-ZGNR/graphene, and the corresponding (b) graphene- 
and (c) nanoribbon-weighted bands, (d)-(f) Corresponding bands for 32-ZGNR/graphene. Bands for the freestanding AFM 
ZGNR are overlaid as blue solid lines for comparison in (b), (c), (e), and (f). Letters L and R stand for the left side and the 
right side of the nanoribbon, respectively. 



FIG. 4. Electronic bands around K. (a)fc-projected bands along F-K (cut A in (e)) and (b) perpendicular to F-K (cut B in (e)) 
for 8-ZGNR/graphene. (c) and (d) corresponding fc-projected bands for 32-ZGNR/graphene. (e) BZs of 8-ZGNR/graphene 
and 1x1 graphene. High symmetry points in each cell are also shown. A and B indicate different cuts about K. (f) schematic 
illustration of half-metallic Dirac cone in the BZ of graphene. 


about K’, which give rise to the same band structures as 
required by symmetry and are not shown here.) Fig. 
schematically show the half-metallic Dirac cones at K 
and K’ in the Brillouin zone (BZ) of graphene. 

A question arises whether the half-metallic behavior 
seen near K is an artifact of the choice of supercell, in par¬ 
ticular the ratio of the graphene substrate to the nanorib¬ 
bon. Although a calculation for an AFM-ZGNR on infi¬ 
nite graphene is obviously infeasible, we performed a cal¬ 
culation for 8-AFM-ZGNR on graphene by doubling the 
ratio, i.e., the ratio of the area of the graphene to that 


of the nanoribbon is increased to 3.5 from 1.75. Such 
a structure gives rise to ~45 Aseparations between the 
nanoribbon and its periodic images on the graphene sub¬ 
strate. Figure shows that there are only minor changes 
in the band structure near K when the ratio is doubled; 
in particular, the half-metallic behavior is well preserved. 
The gap for spin down, however, is slightly reduced, as 
the substrate bands approach those of a single graphene 
layer in the limit of infinite graphene. Nevertheless, one 
may still expect half-metallic behavior as long as there 
exists the spin- and pseudospin symmetry breaking which 
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FIG. 5. Electronic bands along F-X near the point where 
the half-metallic behavior occurs for an 8-AFM-ZGNR on 
graphene with different supercells corresponding the ratio 
(1.75 and 3.5) of the the area of the graphene snbstrate to 
that of the nanoribbon. 

are caused by the presence of the AFM ZGNR and the 
AB stacking. 

Another possible issue is the effect of vdW dispersion 
forces between the nanoribbon and the substrate. To 
investigate such effects vdW-DF calculations were per¬ 
formed for 8-AFM-ZGNR/graphene using the method 
developed by Klimes and Michaelidesp2ESl for which the 
optB88-vdW functional was used for the exchange func¬ 
tional. The equilibrium layer distance derived from these 
calculation is about 3.46 A, close to the 3.35 Aused for the 
above calculations, and thus the calculated band struc¬ 
ture (not shown), is very similar to that shown in Fig. 
including half-metallic nature. 


B. Tight-binding calculations 

The half-metallic states near K originate from the spin 
and pseudospin-symmetry breaking due to the special 
stacking of the two constituents, which can be under¬ 
stood by a tight-binding model. The Hamiltonian for 
the graphene can be written as 

+ h.c.) + (1) 

ij i 

The edge states of ZGNRs can be described by the Hub- 
bard model within the Hartree-Fock approximatiorP^^ 

^ “I" h.C.) 

~\~U ^ ^ 2 ^ ^icr — a ( 2 ) 

i 

i 

where and Cia are creation and annihilation operators 
for spin cr at site i, respectively, Uia = cl^Cia, and tij and 




FIG. 6. TB calculations of AFM-ZGNRs/graphene for 
(a) 8-ZGNR and (b) 32-ZGNR. (c) Band structure for 32- 
ZGNR/graphene when the on-site energies for the top and 
hollow sites are treated differently, and (d) when weak mag¬ 
netism is induced in the substrate. 


U denote hopping integrals and on-site Coulomb inter¬ 
action between electrons, respectively. The respective 
chemical potentials are The self-consistent solution 
to Eq. (§ gives rise to 

^ ^Aa ^ ^B—{T ^ (3) 

with respect to the symmetry center, where A and B are 
the two sublattices corresponding to the different edges. 
Note that Eq. ([^ implies spin degeneracy in freestanding 
AFM ZGNRs. In our calculations U was set to 2.0 eV 
and the tij include up to the second-nearest-neighbor in¬ 
teraction according to Ref [3S] , such that electronic bands 
derived from TB calculations are in good agreement with 
ab initio results for the free standing nanoribbons. 

The Hamiltonian describing the interaction between 
the ZGNR and graphene has the form 

+ h.C.), (4) 

b' 


where a\ and Cj are creation and annihilation operators 
for the two constituents, respectively. Tij describes the 
hopping between atom i in the nanoribbon and atom j in 
graphene. In AB stacking, top-site carbon atoms expe¬ 
rience stronger perturbation than those in hollow sites. 
This physics can be captured by considering hopping only 
between the nearest top-top sites, properly chosen so that 
TB calculations reproduce the bands of graphene bilayer 
about the K point. As a consequence of breaking 
the spin symmetry in Eq.(^, the two edge states of the 
nanoribbon, which carry different spins, are no long de¬ 
generate. Similarly, the interaction breaks the spin and 
pseudospin-symmetry in graphene in the same manner. 

Figures [6}r, b display the electronic bands derived 
from the TB calculations for 8-ZGNR/graphene and 32- 
ZGNR/graphene, respectively. These calculations show 
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FIG. 7. Interaction between bands of AFM-ZGNR and 
graphene, (a) DFT-derived band structures for AFM 32- 
ZGNR/graphene with respect to different interlayer separa¬ 
tions (4.0-5.0 A). Gap openings at k\ and k 2 are labeled by Ai 
and A 2 , respectively, (b) TB bands for the isolated systems 
before interaction. Bands 1 and 2 are the graphene substrate 
linear dispersions. Bands 3 and 4 are the valence and conduc¬ 
tion bands, respectively, of AFM 32-ZGNR. (c) and (d) 714 
and 723 derived from Eq. [^for (c) spin up and (d) spin down. 
The two blue solid lines show the difference between 714 at 
fci and 723 at ^ 2 . 


the same trend as DFT calculations: the gap opening 
for the the majority spin of the top sites is much larger 
than the one for the other spin. If the on-site energies 
for the top and hollow sites are treated differently, as 
commonly done for graphene bilayers, a band shifting at 
the Brillouin zone (BZ) boundaries (the X point) is ob¬ 
tained (Fig. [§:). Because of magnetic proximity effect, 
the coupling of an AFM ZGNR to graphene leads to spin 
polarization in graphene. In this sense, introducing weak 
magnetism in graphene shrinks the gap opening (Fig. §!)■ 
In particular, the gap for spin down becomes extremely 
small, consistent with the DFT calculations. 


C. Low energy model calculations 

To see how the nanoribbon-graphene interaction af¬ 
fects the gap openings at K, we vary the layer distance 
gradually from 5 to 4 A. Because the half-metallic fea¬ 
ture is more pronounced for larger ZGRNs, calculations 
were carried out for 32-ZGNR/graphene only, with the 
results are summarized in Fig. The gap openings at 
K are accompanied by two band splittings at ki (for the 
valence band) and ^2 (for the conduction band) labeled 
as Ai and A 2 , respectively, resulting from the interaction 
between the nanoribbon bands 3 and 4 of Fig. 13^ and the 
linear dispersing bands of the graphene substrate. Based 
on the bands of the isolated systems shown in Fig. [^, 
the band splittings Ai is due to the interaction between 
bands 2 and 3, whereas A 2 is attributed to the interac¬ 
tion between bands 1 and 4. For spin up, which has a 


VB, spin up 

8 ’§'§’& 8 > 8 >&§> 8 >§’ 8 - 8 * 8 « 8 * 8 " 8 » 8 * §• «• s *«»"»“8 "8 ”8 »8 -8 *8 

VB, spin down 

8» & g» 8» 8* 8» 8" 8» 8* «• it ■>» "8 •! >8 «8 •§ .g .g « .g ,g.g ,g<g.g^,g,g 

CB, spin up 

g. g. g. g. g. g. g. 8. §. 8* M M 08 "I *8 >8 «g *8 »g -g <8 € <8 ■€ <8 € ■§ € ■S 

CB, spin down 

§■ S’ & & S’ & 8" 8» 8" 8" 8" §“ 8- 9" 9" 8» M 88 •« »8 "8 -g -B "8 »8 ”8 

FIG. 8 . Band-decomposed charge density at k 2 for an iso¬ 
lated AFM 32-ZGNR. VB and GB denote the valence band 
and the conduction band (labeled by 3 and 4, respectively, in 
Fig.inb). The two sublattices are colored differently. 


sizable gap at K, Ai is much larger than A 2 . For spin 
down, however, A 2 grows and becomes larger than Ai 
as the interlayer separation decreases. Such an enhanced 
splitting pushes the bonding state (resulting from the in¬ 
teraction between bands 1 and 4) down further, reducing 
the gap with the antibonding state resulting from the 
interaction of bands 2 and 3. 

The behavior of Ai and A 2 is inherently related to the 
characteristics of the wave functions of the edge states 
in AFM ZGNR. Unlike in graphene where the two sub¬ 
lattices make equal contributions to the linear dispersing 
bands, they make asymmetric contributions to the va¬ 
lence and conduction bands of AFM ZGNRs. At the 
zone boundary (X point), the valence and conduction 
band edge states are completely localized on only one 
sublattice, either top sites or hollow sites. As the k point 
varies along X-K, these states involve an increasing con¬ 
tribution of the other sublattice but a large asymmetry in 
the contributions of the two sublattices remains. For in¬ 
stance, Fig. 1^ shows that the valence band for freestand¬ 
ing 32-ZGNR near the K point is mainly contributed by 
top sites for spin up and by hollow sites for spin down, 
and for the conduction band the situation is opposite. 
Therefore, one may expect a larger Ai for spin up than 
spin down, as well as a larger A 2 for spin down when the 
coupling to graphene in the AB stacking shown in Fig. 

To gain further insights into the interaction, a 4x4 low 
energy Hamiltonian for the four bands is constructed. 


/ ei 0 Hisik) Hi4{k)\ 

0 £2 H23{k) H24{k) 

HUk) H*^{k) £3 0 

\Ht4{k) m 4 {k) 0 £4 / 


( 5 ) 


where Hij(k) describe the interaction between the bands 
of graphene and those of AFM-ZGNRs shown in Fig. |§5. 
The matrix elements Hij{k) are calculated based on the 
eigenvectors of the freestanding systems: 

H,,{k) =< (fc)|ff“‘|V-f (fc) >= ( 6 ) 

where 'ip'^ and are derived from Eq. Q and Eq. ([^, 
respectively, 7 and 9 stand for the coupling strength and 
phase factor, respectively. 

The band splittings Ai and A 2 are mainly the conse¬ 
quence of H 23 and Hi 4 , respectively. Figures and d 
show the magnitudes of H 14 and H 23 for the two spins, re¬ 
spectively, i.e., 7^4 and 7 ^ 3 . Opposite trends are observed 
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for the two couplings by a comparison of the two plots: 
723 is larger than 723 , but 7 I 4 has a smaller value than 
744 . Thus, one may have a larger Ai for spin up than spin 
dowm as well as a larger A 2 for spin down than spin up. 
Fig. ITc shows that 723 (^ 1 ) is larger than 744 (^ 2 ), while 
Fig.|^ indicates an opposite trend. Moreover, for spin up 
the difference between 723 (^ 1 ) and 714 (^ 2 ) is much larger 
than that for spin down, leading to a noticeable difference 
between Ai and A 2 for spin up, consistent with the ob¬ 
servation in Fig. [7l:a). The 4x4 Hamiltonian reproduces 
the TB band structure shown in Fig. 

In summary, half-metallic Dirac cones are found in AB 
stacking of AFM-ZGNRs/graphene by combining tight- 
binding and DFT calculations. This novel behavior re¬ 
sults from spin and pseudospin symmetry breaking in¬ 
teractions caused by the particular - but common - AB 
stacking. The present finding demonstrates that the 
unique combination of spin and pseudospin in zigzag- 
graphene-nanoribbons can be used to manipulate elec¬ 
tronic properties of graphene. Our results have impli¬ 


cations for both fundamental investigations and practi¬ 
cal applications; like the rich physics related to Dirac 
electrons in graphene, the half-metallic Dirac electrons 
in AFM-ZGNR/graphene may also give rise to extraor¬ 
dinary properties and novel phenomena, such as inter¬ 
esting quantum electronic transport found at graphene- 
monolayer/bilayer junctions.^ Similarly, one may expect 
novel electronic transport properties, with potential ap¬ 
plications in spintronics. We anticipate that our results 
will stimulate further investigations on electronic and 
physical properties of half-metallic Dirac cone in this 
family. 
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